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Using Soft Collinear Effective Theory we analyze the supersymmetric contributions to the branch- 
ing ratios of B s — > 4>ir° and B B — > <f>p° decays. Using mass insertion approximation we show that 
SUSY contributions mediated by chargino exchange can enhance the branching ratios with respect 
to the SM prediction. We show that in single mass insertion scenario, Br (B s —¥ 4>tt°) can be en- 
hanced by about 140% with respect to the SM prediction. In two mass insertion scenario we find 
that Br (B s -> 0tt o ) can be enhanced by about 400% with respect to the SM prediction. For Br 
(B a -t 4>p°) we find that its SM prediction within SCEt is 4.6 x 10" 8 . Including SUSY contribution, 
we show that in single mass insertion scenario Br (B B — > 4>p°) can be enhanced by about 130% with 
respect to the SM prediction. In two mass insertion scenario we find that Br (B s — > (/>7T°) can be 
enhanced by about 160% with respect to the SM prediction. 

PACS numbers: 13.25.Hw,12.60.Jv,11.30.Hv 



I. INTRODUCTION 

The decay modes B —> Ktt, B s — > 07r° and B s — > cj)p are generated at the quark level via b — > s 
transition. Their amplitudes receive contributions from isospin violating electroweak (EW) penguin 
amplitudes. However, these contributions are expected to be small in the case of B — > Ktt that receive 
large contributions from isospin conserving QCD penguins amplitudes which are absent in B s — > (f>TT° 
and B s — > (f>p° decays. Within SM, EW penguin amplitudes are small and hence the predicted branching 
ratios (Br) of B s — > (fin and B s — > <f>p° decays are so small. As a consequence, sizeable enhancement of 
these branching ratios will be attributed only to isospin-violating new physics which can shed light on 
the Ktt puzzle Q] . 

Supersymmetry (SUSY) is one of the best candidates for physics beyond SM. SUSY provides solution 
to the hierarchy problem. Moreover, SUSY provides new weak CP violating phases which can account 
for the baryon number asymmetry and other CP violating phenomena in B and K meson decays [3H5] . In 
addition, the effect of these phases has been studied in the CP asymmetries of r decays in Refs.( [IfjT - H2] ). 

The decay modes B s — > 4>tt° has been studied within SM in framework of QCD factorization in 
Ref.([T3]), in PQCD in Ref.(QH) and in Soft Collinear Effective Theory (SCET) in Ref.([15 ). On the 
other hand, the decay modes B s — > 4>tt° and B s — > (f>p° has been studied within supersymmetry (SUSY) 
in framework of QCDF [21 US]- In this paper we study SUSY contributions to the branching ratios of the 
B s ->■ <pTT° and B s -> <pp° decays in the framework of SCETpIrEO]. 

SCET is an effective field theory describing the dynamics of highly energetic particles moving close 
to the light-cone interacting with a background field of soft quanta[H]. It provides a systematic and 
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rigorous way to deal with the decays of the heavy hadrons that involve different energy scales. Moreover, 
the power counting in SCET helps to reduce the complexity of the calculations and the factorization 
formula provided by SCET is perturbative to all powers in a s expansion. 

In SCET, we start by defining a small parameter A as the ratio of the smallest and the largest energy 
scales in the given process. Accordingly, we scale all fields and momenta in terms of A. Then, the QCD 
lagrangian is matched into the corresponding SCET Lagrangian which is usually written as a series of 
orders of A. The smallness of A allows us to keep terms up to order A 2 in the SCET Lagrangian which in 
turn simplifies the calculations. 

We can classify two different effective theories: SCET/ and SCET// according to the momenta modes 
in the process under consideration. SCET/ is applicable in the processes in which the momenta modes 
are the collinear and the ultra soft as in the inclusive decays of a heavy meson such as B — > A* 7 at the 
end point region and e~p — > e~~ X at the threshold region in which there are only collinear and ultra soft 
momenta modes. SCET// is applicable to the semi-inclusive or exclusive decays of a heavy meson such 
as B — ¥ Dtt, B — ¥ Kir, B —¥ 7we,....etc in which there are only collinear and soft momenta modes. 

This paper is organized as follows. In Sec. [TTJ we briefly review the decay amplitude for B —> MxM 2 
within SCET framework. Accordingly, we analyze the branching ratios of B s — > 4>n° and B s — > cf>p 
decays within SM and SUSY in section |HI| We give our conclusion in Sec. |IV| 

II. B -> AfiM 2 IN SCET 

At leading order in a s expansion, the amplitude of B — > M\Mi where M\ and M 2 are light mesons in 
SCET can be written as 

aSCET _ aLO , AX I Ac.c /in 

/l B->M 1 M 2 — -^B^Mi M 2 + -^B-s-Afi M 2 + •*-B-yM 1 M 2 \ L ) 

Here tf denotes the leading order amplitude in the expansion 1/nib, Ag^ M M2 denotes the 

chirally enhanced penguin amplitude generated by corrections of order a s (/i^)(/iA/A/m 2 ) where /im is 
the chiral scale parameter and •A.b%m 1 m 2 denotes the long distance charm penguin contributions. For 
detail discussion about the formalism of SCET we refer to Refs. [T51 [T7H2U1 12"2"] . 

The decay modes B s — > (j>n° and B s — > (j)p° receive contributions only from Ag°^ MlM2 and so we give 
a brief review for this amplitude in the following. 

At leading power in (1/mt) expansion, the full QCD effective weak Hamiltonian of the = 1 decays 
is matched into the corresponding weak Hamiltonian in SCETj by integrating out the hard scale ntb- 
Then, the SCETj weak Hamiltonian is matched into the weak Hamiltonian SCETjj by integrating out 
the hard collinear modes with p 2 ~ Amf, and the amplitude of the As = 1 decays can be obtained 
via g3]: 
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aLO 

AB->-MiM 2 



-i{MiM 2 \H^ ETlI \B) 

Gf ™ b (j Ml [ [ dudzT Ml j(u, z)Cf M2 {z)4> Ml (u) 



~BM 2 



duT Ml du)4> Ml (u)] + (Mi O M 2 ) 



(2) 



The parameters ^(m 1 ,M 2 ) j (B{M U M 2 ) 

are treated as hadronic parameters and can be determined 
through the % 2 fit method using the non leptonic decay experimental data of the branching fractions and 
CP asymmetries. The hard kernels Ti Ml 
are functions of the Wilson coefficients as follows 22 



Ma )£ and T(M!,M2)J are expressed in terms of and b\ J> which 



df) 



T 1C (u) = Cf L i c[ f) (u) + C™' C% u c 2 /} (u) 
+Cf L M > 4P („) + Cj& c£» («) , 
T U (u, z) = C™* C%\ b[ f) (u, z) + C™> C^ u 4 /} (u, z) 
+Cf L M > bi f \u,z)+ C™> Cj* bf („, z) 



(3) 



here / stands for d or s and C 4 and C,f are Clebsch-Gordan coefficients that depend on the flavor 
content of the final states and and &^ are given by [23] 
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where u; 2 = to^m and W3 = —m^u. u and u = 1 — u are momentum fractions for the quark and antiquark 

( f ) ( f ) 

n collinear fields. The Ac) and Ab\ denote terms depending on a s generated by matching from Hw- 
The 0(a s ) contribution to Ac\ f) has been calculated in Refs.( [2"01 |2"S ] ) and later in Ref. ([22) while 
the C(a s )contribution to Ab[ f) has been calculated in Refs.([22 l l27l l28j ). 



III. NUMERICAL ANALYSIS AND RESULTS 



In this section, we analyze the branching ratios for B 8 — > (fiir and B s — > <j>p° decays. We take 
m& = 4.7 GeV and the Wilson coefficients Ci at leading logarithm order are are given by [2T)j : 

Ci_io(m 6 ) = {1.110 , -0.253 , 0.011 , -0.026 , 0.008 , -0.032 , 0.09 x 10~ 3 , 0.24x 10~ 3 , 

-10.3xl0~ 3 ,2.2xl0~ 3 }. (6) 
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For the other hadronic parameters, we use the same input values given in Ref . ( [2"2"] . 
At quark level, the decay modes B s — > <fiir° and B s — > <fip° are generated via b — > s transition and 
hence we can decompose their amplitudes A according to the unitarity of the CKM matrix as 

A(B S -> </m°,<j>p ) = K{AT e + A% w ) + \ s c A™ n - cc (7) 

Here A^ = V p bV* s with p = u,c and A^[ ee ,A^ w refer to the color suppressed tree and Electroweak 
penguins amplitudes respectively. JS^° n ~ cc refers to contributions from Electroweak penguins which are 
proportional to A*. In the SM, without including QCD corrections, we find that A% ee > A% w , A n c on - cc 
due to the hierarchy of the Wilson coefficients C'2 3> C3— 10 • Thus, we can write, to a good approximation, 

A(B S -> <M°,0p°) ~ \ s u A% ee + \ s c AT n -° c (8) 

which shows as A^ ee will be suppressed by a factor |A£| ~ 0.02 1 A* | in comparison to A^ on ~ cc . Thus we 
see that Electroweak penguins amplitudes can compete with the tree color suppressed amplitude and one 
expects that new physics contribution to the penguin amplitude, A™ on ~ cc , can enhance significantly the 
decay rates and the branching ratios of these decay modes. 

Another remark here is that since these decay modes do not receive contributions from the long distance 
charm penguin and so we expect very small branching ratios as non-perturbative charming penguin plays 
crucial rule in the branching ratios using SCET. 

A. B s -> (j>ir° 

At leading order in a s expansion the amplitude of B s — } <fiir° decay including SUSY contributions is 
given by 

A(B° S -> </>7T°) x 10 6 ~ (-3.6C10 + lTCio + 8.3C 7 - 8.3C 7 + 1.9C 8 - l.9C s - 8.3C 9 + 6.6C , 9 )A ts 

+ (2.4Ci - 0.9Ci + 5.6C2 - 4.4C*2)A US (9) 

Where Ci and Ci are the Wilson coefficients which can be expressed as 

d = Cf M + Cf USY , Ci = Cf USY (10) 

Ci are generated from the weak effective Hamiltonian by flipping the chirality left to right and so in the 
SM we have Cf — 0. The dominant SUSY contribution to the Wilson coefficient come from diagrams 
with exchanging gluino and chargino and so we can write 



Cf USY = Cf + Cl (11) 

where Cf represents the gluino contribution and Cf represents the chargino contribution. 

In SUSY, Flavor Changing Neutral Current(FCNC) and CP quantities are sensitive to particular 
entries in the mass matrices of the scalar fermions. Thus it is useful to adopt a model independent- 
parametrization, the so-called Mass Insertion Approximation (MIA) where all couplings of fermions 
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FIG. 1: Br (_B S — > 4>tt°) x 10 6 versus the phase of the mass insertion (<5^ L )32. The horizontal line in the diagram 
represents the SM prediction. 

and sfermions to neutral gauginos are flavour diagonal [3U]. The complete expressions for the gluino 
and chargino contributions to the Wilson coefficients in MIA can be found in Refs. [3l l31~M33| . After 
substituting with Wilson coefficient and keeping only dominant terms, we get 



A{B S -> (jjit") x 10 y ~ -2.781 - 1.005z - 0.006(c52 L ) 23 + 0.005(^ fl ) 23 + 0.014(££ L ) 32 + 0.334(££ R ) 



31 



+ 1.520(^) 32 - 1.730(5^)32 - 0.381(^)31 - 0.015(^)32 (12) 

Setting all mass insertions to zero which corresponds to SM case, we find Br(B s —> <fiir°) = 7 x 10 -8 
which is in agreement with prediction in Ref.([15 ). 

SUSY contribution to Br(B s — > </>7r°) requires non vanishing mass insertions which in general are 
complex and thus they have new sources for weak phases that can lead to unwanted effects in CP 
violation phenomena. Therefore, these mass insertions are subjected to constrains from vacuum stability 
argument [34] , experimental measurements concerning FCNC and CP violating phenomena [3]. Recent 
studies about other possible constraints can be found in Refs.( 



As can be seen from Eq.(12) the dominant SUSY contributions are those proportional to the mass 
insertions {8^)32 , {^lr)^ which come from chargino contributions to Wilson coefficients. Thus we can 
approximately write 



A(B S 07r ) x 10 9 ~ -2.781 - 1.005i + 1.520(<^) 32 - 1.730(<5£ x ) 32 (13) 

The strongest constraints on the mass insertions (8rl)32 and (££ fl ) 32 are imposed by vacuum stability 
argument [3H |3S]. Applying vacuum stability argument, we find that l^^al = 1^^)32! — C(l) for 
squark masses < 700 GeV and hence to maximize SUSY contributions we can parameterize the mass 
insertions as (§%l)32 = e l5s and (6)^)32 = e l&L where 5r and Sl are the phases of the mass insertions 
( S rl)32 and (dl R ) 32 respectively. 
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FIG. 2: Left diagram corresponds to Br (B s — > 10 versus the phase of the (^£_b)32 where we assume 

| (5^^)32 1 — 0(1) and the horizontal line represents the SM prediction. Right diagram corresponds to contours of 
Br (B s — > 07T°) x 10 6 versus the phases of the (5^)32 and (8^)32 for two mass insertions scenario. 

In our analysis we consider two scenarios, the first one with a single mass insertion where we keep only 
one mass insertion per time and take the other mass insertion to be zero and in the second scenario we 
consider non vanishing two mass insertions. 

We start by considering the first scenario and set (8^)32 = 0. Thus, in this case, the amplitude in 



Eq.( 13 ) becomes 

A(B° S 07r°) x 10 9 ~ -2.781 - 1.005* - 1.730(c5^ L ) 32 (14) 



The phase of the mass insertion (S l ^ L )s2, 8r, can be varied from — tt — > tt. Clearly from Eq.(14|, last 
term will have opposite sign to the other terms for the values 8r — ±tt which leads to a destructive 
contribution to the amplitude and thus we find Br (B s — > 4>tt°) ~ 1.6 x 10~ 8 which is less than the SM 



prediction. For a value <5r = last term in Eq.( 14) will have same sign as other terms. As a consequence, 
the amplitude will be enhanced and we find that Br (B s —> 4>tt°) « 1.7 x 10~ 7 which is more than 100 % 
enhancement of the SM prediction. Variation of the Br (B s — » <fiir°) with 8r is plotted in Fig. ([I]) where 
the horizontal line represents the SM prediction. As can be seen from Fig.([l]) that, Br (B s — > 4>tt°) is 
enhanced by SUSY contribution for Sr lies in the range ] — 1.5, 2. 2 [ where the enhancement can reach 
about 140 % with respect to the SM prediction. 

Turning now to the case {5^ R )32 7^ and set the other mass insertion to zeros. In this case the 
amplitude becomes 

A(B° S 07T ) x 10 9 ~ -2.781 - 1.005* + 1.520(^ R ) 32 (15) 



As can be seen from Eq.( 15 ), last term will have a sign similar to the other terms for the values Sl — ±1" 



which leads to a constructive contribution to the amplitude and thus we find Br (B s — > (fin ) w 1.55 x 10 7 



which is is more than 100 % enhancement of the SM prediction. For a value 8l — 0, last term in Eq.(15l 
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FIG. 3: Br (B a 
SM prediction. 



o°) x 10 s versus the phase of the {8lr)s,2. The horizontal line in the diagrams represents the 



will have opposite sign as the other terms. As a consequence, the amplitude will be smaller than the SM 
case and we find that Br (B s — > </>7r°) 2.1 x 10 -8 which is less than the SM prediction. Variation of the 
Br (B s — !• (f>TT°) with Sl is plotted in Fig.([2| left where, as before, the horizontal line represents the SM 
prediction. As can be seen from Fig.Q left that, Br (B s — > <fiir°) is enhanced by SUSY contribution for 
all values of Sl except in the range [—1, 1.8]. 

We consider now the second scenario by considering two mass insertions per time. As can be seen from 
Eq.(13) for Sl — ±7r and Sr = last two terms will have same sign as the other terms and this will 
enhance the amplitude and thus we find that Br (B s <fiir°) ~3x 10 -7 which is enhancement by about 
400 % with respect to the SM prediction. For Sl = and Sr = ±7r last two terms will have opposite sign 
to the other terms and this will decrease the amplitude and thus we get Br (B s — > 4>tt°) rjIx 10~ 8 . In 
FigQ right, we plot the contours of Br (B s (f>ir°) as a function of 5l and Sr. As can be seen from that 
Figure, Br {B s — > 4>ir°) is enhanced in the regions of the parameter space where the two mass insertions 
have constructive contributions and is decreased in the regions of the parameter space where the two 
mass insertions have destructive contributions. 



B a 



The decay mode B® — > <j>p° contains two vector mesons in the final state and thus it is characterized by 
three helicity amplitudes Aq (longitudinal) and A±. Naive factorization analysis leads to the hierarchy 
A : A- : A + = 1 : A- : (^) 2 where m& ~ 5GeV is the bottom quark mass and A rj 0.5GeV is the strong 
interaction scale [38 a . The hierarchy shows that the dominant contribution is mainly from the longitudinal 
polarization component. Thus in our calculation we consider only the longitudinal amplitude and so we 
get 
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FIG. 4: Contours of Br (B s -t <j>p°) x 10 6 for two mass insertions scenario. Left diagram corresponds to (8^11)32 
and (82r)31- Right diagram corresponds to (5^)32 and (5^)32 



A(B S -> </>p°) x 10 6 ~ (-4.9Cio - 4.9<5i ~ 16.5C 7 + 16.5CV - 2.2C 8 + 2.2C* 8 - 4.6C 9 - 4.5C 9 )A ts 

+ (3.3C X + 3.2<5i + 3.1C 2 + 3.0<5 2 )A„ S (16) 

We obtain after substituting with Wilson coefficients 



A(B S -> ^p ) x 10 9 



-0.445 - 2.353 4 - 0.003(^ L ) 23 - 0.003(4 R ) 23 - 0.031(<5£ L ) 32 - 0.293(J^) 3 i 



1.330(5^)32 + 0.113(^)32 + 0.025(^)31 + 0.030(^ fl )32 



(17) 



Setting all mass insertions to zero which corresponds to SM case, we find Br(B s — > 4>P°) = 4-6 x 10 8 . 
As can be seen from eq.(fl7l) the dominant SUSY contributions are those proportional to the mass 



insertions (6^)32 , (^2^)32 and (5^)31 ■ Thus we can rewrite the amplitude in Eq.(17) as 

A(B S ->■ <f>p°) x 10 9 ~ -0.445 - 2.353i - 0.293(<J£ fl ) 3 i - 1.330(^ H ) 3a + 0.113(<J£ L ) 3 2 (18) 



We repeat the same analysis as in subsection (III A) by considering two different scenarios one with 



single mass insertion per time and the other one with two mass insertion per time. We start by considering 
the first scenario and set all mass insertions to zero except the mass insertion (#£r)s2- Thus, in this case, 



the amplitude in Eq.(18) becomes 
A(B S -> 4> 



>°) x 10 9 



-0.445 - 2.353i - 1.330(<^) 3 2 



(19) 



As can be seen from Eq.(19), last term will have a sign similar to the other terms for the value Sl = 
which leads to a constructive contribution to the amplitude and thus we find Br (B s — > <fip°) ~ 7 x 10~ 8 
which shows that SUSY contribution enhances the SM prediction by about 50%. For the values Sl = ±tt, 



last term in Eq.(19l will have opposite sign to the other terms. As a consequence, the amplitude will 
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be smaller than the SM case and we find that Br (B s — > (ftp ) ~ 5 x 10~ 8 which is 9% enhancement to 
the SM prediction. Variation of the Br (B s — > <j>p°) with Sl is plotted in Fig.([3| where, as before, the 
horizontal line represents the SM prediction. As can be seen from the Figure, Br (B s — > (ftp ) can be 
enhanced by about 130% with respect to the SM prediction. 

Turning now to the case (S RL )32 7^ and other mass insertions equal zero. We find in this case, 
Br (B s — > (ftp ) can be enhanced only by 5%. We consider now the case (6^)31 ^ and other mass 
insertions equal zero. It should be noted that single-top production involves the mass insertion (<5]^)3i 
and thus direct top production is the only way to observe or constrain (<5"^)qi |55J. With unconstrained 
(^Ei?j3i we nn d that the enhancement can reach 20%. 

We consider now the second scenario by considering two mass insertions per time and set other mass 
insertion to zero. We consider first the two mass insertions (Sf /R )32 and (S l £ R )3i and set {8 RL )32 — 0. In 



this case the amplitude in Eq.(18) is reduced to 



A{B S (ftp ) x 10 9 ~ -0.445 - 2.353i - 0.293(^)31 - 1.330(££ R ) 32 ( 20) 



Clearly from Eq.(20), Br (B s — » (ftp ) will be enhanced significantly when both (^ H ) 3 2 and (Sl R )3i 
have phases > which is shown in Fig.Q left where we plot the contours of Br (B s — > (ftp ) versus the 
phases of the (S^ R )32 and (S^r)^- ^ n addition, Br (B s — > (ftp ) becomes less than SM prediction when 
both (<*>2fl)32 and (S^ R )3i have negative phases. Finally, we see also from the Figure, two mass insertions 
enhance the branching ratio by about 160% with respect to the SM prediction which is larger than the 
enhancement obtained in the first scenario with one mass insertion. 

Finally, we plot in Fig.Q right Br (B s — > (ftp ) as a function of the phases of the two mass insertions 
{81 R )32 and {5 RL )32 where we set {5 RL )3i = 0. As can be seen from the Figure, Br {B s — > (ftp ) is also 
enhanced by about 100% with respect to the SM prediction although the enhancement in this case is 
smaller than the previous case with the two mass insertions (Sf iR )32 and (5^ R )3i- 



IV. CONCLUSION 



We have analyzed SUSY contributions to the branching ratios of (B s — > (ftir ) and (B s — > (ftp ) decays 
using Soft Collincar Effective Theory. We have adopted in our analysis the mass insertion approximation 
in two different scenarios. In the first scenario we keep one mass insertion per time and set others to 
zero while in the second scenario we keep two mass insertions and set others to zero. For squark masses 
< 700 GeV , where the constraints on the mass insertions in the up sector are not effective, we have 
shown that, Br (B s — > (ftir ) can be enhanced by about 140% with respect to the SM prediction in the 
first scenario. In two mass insertion scenario we find that Br (B s — > (ftir ) can be enhanced by about 
400% with respect to the SM prediction. 

For Br (B s (ftp ) we find that its SM prediction within SCEt is 4.6 x 10~ 8 . Including SUSY 
contribution, we have shown that Br {B s — > (ftp ) can be enhanced by about 130% with respect to the 
SM prediction in the first scenario. In the second scenario we find that Br (B s —> (ftir ) can be enhanced 
by about 160% with respect to the SM prediction. 
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